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Abstract. We present Sunyaev-Zel'dovich (SZ) effect observations of a sample of 25 
massive relaxed galaxy clusters observed with the Sunyaev-Zel'dovich Array (SZA), an 
8-element interferometer that is part of the Combined Array for Research in Millimeter- 
wave Astronomy (CARMA). We perform an analysis of new SZA data and archival 
Chandra observations of this sample to investigate the integrated pressure — a proxy 
for cluster mass — determined from X-ray and SZ observations, two independent probes 
of the intra-cluster medium. This analysis makes use of a model for the intra-cluster 
medium introduced by Bulbul (2010) which can be applied simultaneously to SZ and 
X-ray data. With this model, we estimate the pressure profile for each cluster using a 
joint analysis of the SZ and X-ray data, and using the SZ data alone. We find that the 
integrated pressures measured from X-ray and SZ data are consistent. This conclusion 
is in agreement with recent results obtained using WMAP and Planck data, confirming 
that SZ and X-ray observations of massive clusters detect the same amount of thermal 
pressure from the intra-cluster medium. To test for possible biases introduced by our 
choice of model, we also fit the SZ data using the universal pressure profile proposed 
by Arnaud (2010), and find consistency between the two models out to rsoo in the 
pressure profiles and integrated pressures. 
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1. Introduction 



The Sunyaev-Zel'dovich (SZ) effect ( jSunyaev fc Zerdovich , 1972) is a spectral distortion 
of the cosmic microwave background (CMB) caused by the scattering of CMB photons 
off the hot electrons of the intra-cluster medium (ICM). Over the past two decades, SZ 
observations with both single-dish and interferometric instruments have become routine 



(e.g., Birkinshaw et al., 1991; Carlstrom et al., 1996; Holzapfel et al., 1997; Carlstrom 



et al., 2002), and SZ surveys are now producing catalogs of newly-discovered clusters 



out to high redshift (Vanderlinde et al., 2010; Marriage et al., 2011; Williamson et al. 



2011; Planck Collaboration et al., 2011a). SZ measurements are complementary to the 



X-ray measurements which have long been used to study clusters, but only in recent 
years have sufficiently large samples of objects been observed in the SZ to permit a 



rigorous comparison between these two techniques (e.g., Reese et al., 2002; Bonamente 



et al. , 2006; LaRoque et al. , 2006). 



The SZ effect causes a perturbation AT of the CMB temperature Tqmb given by 



AT 



CMB 



M = f{x)y, 



(1) 



where f{x) is the frequency dependence of the SZ effect (e.g., LaRoque et al., 2006); ay 
is the Thomson cross-section; ng, Tg and rrie are the number density, temperature, and 
mass of the electrons, respectively; k is the Boltzmann constant; c the speed of light; 
and the integral is along the line of sight £. At a given frequency, the amplitude of the 
effect depends linearly upon the Compton y parameter, which is defined implicitly in 
Equation [T] Note that the y parameter is proportional to the ICM pressure integrated 
along the line of sight. At frequencies below 218 GHz, the SZ effect causes a decrement 
in the CMB temperature in the direction of the cluster. The integral of y over the 
solid angle Q subtended by the cluster, known as the (cylindrical) integrated Compton 
y parameter l^cyi = Jy dQ^ is expected to be a good proxy for cluster total mass since 
it traces the thermal energy content of the cluster gas. Alternatively, the Compton y 
parameter can be integrated spherically, 

1 f kaT\ 
D\ \meC^ J 

where the volume 1/ is a sphere centered on the cluster and Da is the angular diameter 
distance. 

X-ray data can also be used to constrain the density and temperature — and thus 
the pressure — of the ICM. Over the past decade, several groups have investigated the 
consistency between X-ray and SZ pressure measurements. Early measurements of the 



^sph(^50o) 



ripTpdV 



(2) 



SZ signal from WMAP by, e.g.. Lieu et al. (2006), Bielby & Shanks (2007), detected an 



SZ signal at a lower level than expected. Atrio-Barandela et al. (2008) showed that the 



isothermal beta model leads to an electron pressure profile that exceeds the measured 
values at large radii by a factor of few, and that the baryon profile is consistent with a 



model based on the Navarro et al. ( 1997| ) matter profile. Diego & Partridge (2010) also 
showed that contamination by compact radio sources may have led to underestimates 
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of the SZ effect flux decrements in the WMAP data. More recent comparisons of 



Chandra X-ray data to stacked data from WMAP and Planck (Mehn et al. , 2011 ; Planck 



CoUaboration et al. , 2011b) have found consistency between SZ and X-ray measurements 



for large samples of clusters. Komatsu et al. (2011) also analyzed a sample of massive 
nearby clusters individually resolved by WMAP^ again flnding good agreement with 
X-ray predictions. 



In this paper, we present Sunyaev-Zel'dovich Array (SZA) observations of the Allen 



et al. (2008) sample of massive relaxed galaxy clusters, together with archival Chandra 



X-ray observations that are available for all clusters in this sample. The sensitivity 
and resolution of our data permit us to measure the pressure proflle and the integrated 
pressure out to rsoo — the radius within which the average cluster density is 500 times 
the critical density — for each cluster individually, without the need to resort to scaling 



relations between the X-ray luminosity and mass (as was done by Melin et al., 2011 



Planck Collaboration et al., 2011b, for example). We use the Bulbul et al. (2010) model 



of the cluster pressure, density, and temperature. Since this model has a consistent 
parameterization for all thermodynamic quantities, it is especially well-suited for joint 
X-ray and SZ analysis. As a cross-check against model-dependent biases, we also flt the 
SZ data using the model of jArnaud et al. (2010) based on the numerical simulations 
of Nagai et al. (2007). We flnd consistency to within our measurement uncertainties 
both between the X-ray and SZ measurements, and between the Bulbul et al. (2010) 
and Arnaud et al. (2010) models. 

The paper is structured as follows: Section [2] describes our observations and our 
sample. Section |3] presents our joint analysis technique. Section [4] describes our method 
of measuring the integrated l^sph(^5oo) parameter (deflned in Equation [2]), Section [5] 
presents and discusses our results, and our conclusions are presented in Section |6| 



2. Observations 



The SZA is an eight-element interferometer designed to detect and image the SZ effect 
from clusters at z> 0.1, and is part of the Combined Array for Research in Millimeter- 
wave Astronomy (CARMA). The array is equipped with 30 and 90 GHz receivers; all 
SZA observations presented in this paper were taken at 30 GHz. At this frequency, the 
3.5 m diameter SZA telescopes have a fleld-of-view (or primary beam) of 10. ''T FWHM. 
Interferometric data are proportional to the Fourier transform of the sky brightness. 
These visibility data, denoted V{u^v)^ are sampled at Fourier-plane coordinates (u^v) 
corresponding to the projected separation of pairs of telescopes (or baselines)^ as viewed 
by the source at the time of observation. At the time of the observations discussed in 
this work, the SZA antennas were arranged in a hybrid conflguration, with six closely 
spaced telescopes and two "outriggers" located ^50 m from the inner array. The inner 
six telescopes probe small {u^ v) Fourier modes, sampling the angular scales where the 
SZ signal is largest for moderate- to high-redshift clusters (1 — 6'). Baselines involving 
the outriggers are sensitive to angular scales down to ^ 20^' and are used to constrain 
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the positions and fluxes of unresolved radio sources. 

Of the 42 clusters in the [Allen et al. (2008) sample of massive relaxed galaxy 
clusters, the SZA has observed the 31 objects above 5 > —15° at redshift z > 0.09. 
The declination restriction is imposed by the latitude of the observatory in the Owens 
Valley, California (37°14'02''N, 118°16'56'W), while the exclusion of clusters at low 
redshift is due to the inability of an interferometer to constrain scales larger than 
that which the shortest antenna spacing can probe at the lowest frequency band. The 
largest angular wavelength measured by the SZA is 10.9', which for massive low-redshift 
clusters is generally smaller than 2r5oo/^A- Of these 31 clusters observed with the 
SZA, Abell 2390 and Abell 611 were excluded from this analysis because they did not 
have available local background in their Chandra ACIS-S X-ray observations. Three 
additional clusters — 3C295, C1J1415. 2+3612, and Abell 963 — were discarded because 
of extended or otherwise diflicult-to-remove radio source contamination, and one — 
RXJ0439. 0+0521 — because of a pointing error. 

Our sample therefore consists of 25 clusters. The synthesized beam of the long 
(short) baseline data for this sample is approximately 15-30^' (90-180^''), and the average 
rms noise in the maps is ^ 0.25 — 0.30 mJy. In all cases, the Chandra data provide 
spatially resolved X-ray spectroscopy and sub-arcsecond imaging. A summary of the 
data is provided in Table [T] 

Radio sources detected in the cluster flelds are reported in Table [2j For each cluster 
field, we use the NRAO VLA Sky Survey (NVSS) and Faint Images of the Radio Sky 
at Twenty- centimeters (FIRST) 1.4 GHz catalogs as a reference for locating compact 
radio sources within 10' of the cluster center. Most radio sources in our observations 
have counterparts in the FIRST survey, which has an rms noise of 0.15 mJy at 1.4 GHz. 
Inverted spectrum sources that may be present at 30 GHz may not have counterparts 
at 1.4 GHz, but fortunately they comprise a small fraction of the source population 
dMuchovej et al.[|2010D . 

For all 25 clusters in our sample we have available archival Chandra X-ray 
observations (Allen et al., 2008). Event files for all cluster observations and additional 
blank-sky composite event files used for background subtraction were reduced using 
CIAO 4.3.1 and CALDB 4.3. X-ray spectra are extracted in several annular regions 
for each cluster, centered at the peak of the X-ray emission. Emphasis is placed on 
the removal of periods of high background, and on the modeling of soft X-ray residuals 
that may be present after the subtraction of the blank-sky background. The method of 
analysis of the Chandra data and examples of the temperature and surface brightness 
profiles can be found in Bulbul et al. (2010) and Hasler et al. (2011). More details on 
the Chandra data for all clusters in this sample will be shown in a forthcoming paper, 
in which we will present the measurement of the gas mass fraction from the X-ray 
observations (Hasler et al. in prep.). 

In Figure [l} we show the raw Chandra X-ray images (binned in the 0.7-7 keV 
energy band) for each of the 25 clusters, with contours obtained from the short baseline 
point source-removed SZA data overlaid. 
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Table 1. Sample of massive and relaxed clusters from the Allen et al. (2008) sample 



with high resolution Sunyaev-Zel'dovich effect SZA observations. 



Cluster 


z 




R.A. 


Dec. 


Nh 


SZA 




Chandra 








(J2000) 


(J2000) 


(10^0 cm-2)« 


(hrs)^ 


ACTS 


ObsID 


(ks)^ 


MACSJ0159.8-0849 


0.40 


01 


59 49.5 


-08 50 02 


2.06 


21.2 


I 


3265 


16.4 
















I 


6106 


34.3 
















I 


9376 


19.5 


Abell 383 


0.19 


02 


48 03.4 


-03 31 44 


3.40 


25.0 


I 


524 


9.9 
















I 


2320 


18.5 


MACSJ0329. 7-0212 


0.45 


03 


29 41.7 


-02 11 48 


3.43 


8.1 


I 


6108 


32.7 
















I 


3257 


9.6 
















I 


3582 


19.3 


Abell 478 


0.09 


04 


13 25.2 


+10 27 52 


34.29 


37.1 


I 


6102 


10.0 


MACSJ0429.6-0253 


0.40 


04 


29 36.1 


-02 53 08 


3.23 


22.1 


I 


3271 


23.2 


3C186 


1.06 


07 


44 17.5 


+37 53 17 


5.11 


13.7 


S 


9407 


66.3 
















s 


9408 


39.6 
















s 


9774 


75.1 
















s 


9775 


15.9 


MACSJ0744.9+3927 


0.69 


07 


44 52.9 


+39 27 26 


5.66 


12.7 


I 


6111 


49.5 
















I 


3197 


20.2 
















I 


3585 


19.7 


MACSJ0947.2+7623 


0.34 


09 


47 13.1 


+76 23 14 


2.28 


11.5 


I 


2202 


11.7 


Zwicky 3146 


0.29 


10 


23 39.6 


+04 11 10 


2.46 


6.8 


I 


909 


45.2 
















I 


9371 


36.3 


MACSJ1115.8+0129 


0.35 


11 


15 52.0 


+01 29 58 


4.34 


26.2 


I 


9375 


39.6 


MS1137.5+6625 


0.78 


11 


40 22.2 


+66 08 14 


0.95 


19.6 


I 


536 


109.6 


Abell 1413 


0.14 


11 


55 18.2 


+23 24 19 


3.60 


43.2 


I 


5003 


66.6 
















I 


1661 


9.1 
















I 


5002 


34.4 


C1J1226. 9+3332 


0.89 


12 


26 58.2 


+33 32 47 


1.83 


16.0 


I 


5014 


31.6 
















I 


3180 


29.9 


MACSJ1311.0-0311 


0.49 


13 


11 01.7 


-03 10 38 


1.82 


4.6 


I 


6110 


63.0 
















I 


3258 


13.1 
















I 


9381 


29.0 


RXJ1347.5-1145 


0.45 


13 


47 30.6 


-11 45 10 


4.60 


25.7 


I 


3592 


54.8 


Abell 1835 


0.25 


14 


01 02.0 


+02 52 40 


2.04 


9.0 


I 


6880 


117.9 
















I 


6881 


36.8 
















I 


7370 


40.0 


MACSJ1423.8+2404 


0.54 


14 


23 47.9 


+24 04 42 


2.20 


6.5 


I 


1657 


18.2 


MACSJ1427.3+4408 


0.49 


14 


27 16.3 


+44 07 29 


1.19 


17.4 


I 


6112 


8.8 
















I 


9380 


25.8 
















I 


9808 


14.9 


RXJ1504. 1-0248 


0.21 


15 


04 07.5 


-02 48 16 


5.97 


9.2 


I 


5793 


39.2 
















I 


4935 


11.9 


MACSJ1532.9+3021 


0.36 


15 


32 53.8 


+30 20 58 


2.30 


14.6 




1665 


8.2 


MACSJ 1621.6+3810 


0.46 


16 


21 24.9 


+38 10 08 


1.13 


44.0 




6172 


29.2 


















3254 


9.6 


















6109 


36.7 


















9379 


29.7 


















10785 


29.7 


Abell 2204 


0.15 


16 


32 46.9 


+05 34 31 


5.67 


19.6 




7940 


76.9 


MACSJ1720.3+3536 


0.39 


17 


20 16.8 


+35 36 25 


3.46 


36.2 




6107 


29.1 


















3280 


20.6 


















7718 


7.0 


RXJ2129.6+0005 


0.23 


21 


29 40.0 


+00 05 18 


3.63 


24.5 




552 


10.0 


Abell 2537 


0.29 


23 


08 22.2 


-02 11 28 


4.62 


24.8 




9372 


38.5 








a: Nh 


is HI Galactic column density. 











b: SZA exposure is unflagged, on-source time 
: Chandra exposure is unflagged, on-source time 
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Table 2. SZ Centroids and Radio Source Locations for the SZA Observations 







SZ Centroid 




30 GHz Source 




1.4 GHz Fl 


ux (mJy) 


Cluster 


z 


a(J2000) 


(5(J2000) 


src 






Flux (mJy) 


NVSS 


FIRST 


MACSJ0159.8-0849 


0.40 


01:59:51.5 


-08:50:06.9 


1 


-31.3 


8.0 


84.4+0.2 


36.7 


31.4 


Abell 383 


0.19 


02:48:03.5 


-03:31:55.8 


1 


-1.5 


11.5 


4.3+0.2 


40.9 


- 










2 


276.7 


-149.2 


7.5+0.3 


54.9 


_ 


MACSJ0329.7-0212 


0.45 


03:29:40.3 


-02:11:44.5 


1 


263.2 


-97.2 


12.7+0.4 


37.2 


- 


Abell 478 


0.09 


04:13:25.0 


+10:27:50.8 


1 


195.4 


18.1 


2.9+0.1 


47.7 


_ 










2 


3.3 


4.1 


2.3+0.1 


36.9 


_ 


MACSJ0429.6-0253 


0.40 


04:29:35.6 


-02:53:01.6 


1 


7.0 


-7.7 


18.2+0.2 


138.8 


_ 










2 


-285.4 


53.9 


3.2+0.2 


18.3 


_ 


3C186 


1.06 


07:44:14.8 


+37:53:21.2 


1 


40.8 


-3.3 


22.6+0.2 


1236.4 


1244.9 










2 


81.6 


-94.3 


10.6+0.2 


105.4 


49.2 


MACSJ0744.9+3927 


0.69 


07:44:52.2 


+39:27:34.6 


1 


-215.6 


286.7 


2.5+0.5 


- 


4.4 


MACSJ0947.2+7623 


0.34 


09:47:12.4 


+76:23:03.0 


1 


11.1 


11.2 


2.5+0.3 


21.7 


_ 


Zwicky 3146 


0.29 


10:23:38.9 


+04:11:27.7 


1 


92.4 


-48.5 


5.0+0.3 


95.8 


56.7 










2 


8.9 


-17.6 


1.7+0.2 


7.1 


2.0 










3 


-46.2 


-138.7 


2.2+0.2 


31.5 


15.1 


MACSJ1115.8+0129 


0.35 


11:15:52.2 


+01:29:50.6 


1 


128.7 


-360.5 


2.7+0.2 


11.5 


10.5 










2 


162.2 


-95.5 


3.7+0.2 


_ 


2.1 










3 


139.7 


-64.8 


1.8+0.3 


- 


- 










4 


-3.1 


3.0 


1.4+0.4 


6.2 


5.6 










5 


-168.9 


-71.2 


2.7+0.2 


10.3 


6.4 


MS1137.5+6625 


0.78 


11:40:22.8 


+66:08:13.2 


_ 


_ 


_ 


_ 


_ 


_ 


Abell 1413 


0.14 


11:55:17.5 


+23:24:04.0 


1 


-117.0 


135.5 


2.1+0.1 


28.1 


19.8 










2 


-386.0 


-185.2 


2.8+0.3 


_ 


_ 


CLJ1226.9+3332 


0.89 


12:26:57.7 


+33:32:51.8 


1 


263.4 


-46.0 


3.9+0.2 


29.8 


23.2 


MACSJ1311.0-0311 


0.49 


13:11:02.2 


-03:10:47.0 


- 


_ 


_ 


- 


- 


- 


RXJ1347.5-1145 


0.45 


13:47:31.4 


-11:45:16.1 


1 


-11.4 


6.3 


8.7+0.2 


45.9 


_ 










2 


-339.1 


251.5 


6.9+0.4 


365.8 


- 










3 


-53.3 


279.5 


2.6+0.3 


5.1 


_ 


Abell 1835 


0.25 


14:01:02.2 


+02:52:34.4 


1 


-1.5 


9.2 


2.9+0.3 


39.3 


31.3 










2 


-29.0 


-47.7 


1.0+0.3 


_ 


_ 


MACSJ1423. 8+2404 


0.54 


14:23:48.6 


+24:05:13.6 


1 


-11.8 


-31.5 


2.0+0.2 


8.0 


5.2 


MACSJ1427.3+4408 


0.49 


14:27:15.8 


+44:07:41.4 


1 


4.8 


-10.8 


16.4+0.2 


47.9 


41.3 










2 


33.6 


206.8 


1.1+0.2 


8.6 


8.2 


RXJ1504. 1-0248 


0.21 


15:04:07.1 


-02:48:17.8 


1 


5.7 


1.3 


15.9+0.2 


60.5 


40.8 


MACSJ1532.9+3021 


0.36 


15:32:54.0 


+30:20:59.0 


1 


-39.3 


-72.9 


5.7+0.2 


7.9 


6.0 










2 


-2.7 


0.3 


3.2+0.2 


22.8 


15.2 










3 


-82.8 


-128.4 


1.3+0.2 


18.0 


4.1 


MACSJ1621.6+3810 


0.46 


16:21:25.3 


+38:09:56.9 














Abell 2204 


0.15 


16:32:47.2 


+05:34:34.7 


1 


-3.6 


-1.5 


7.0+0.2 


69.3 


57.9 










2 


-421.8 


-362.8 


21.6+0.2 


41.6 












3 


191.0 


-132.8 


0.7+0.1 


12.2 


1.2 


MACSJ1720.3+3536 


0.39 


17:20:16.2 


+35:36:36.0 


1 


650.3 


340.2 


167.7+0.2 














2 


10.3 


-9.4 


1.8+0.4 


18.0 


16.8 


RXJ2129.6+0005 


0.23 


21:29:40.2 


+00:05:20.9 


1 


-3.2 


0.4 


2.6+0.2 


25.4 


23.8 










2 


228.1 


160.9 


3.1+0.2 


34.3 


6.6 


Abell 2537 


0.29 


23:08:19.2 


-02:11:19.0 


1 


138.6 


437.2 


8.4+0.9 


69.9 


58.6 



a: Offset from fit SZ Centroid 
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2^0"'10' 2^0"'0'' 50' 40' 30' 2V8'"20' 10' 48"'0' 50' 40' S^SO^'O' 50' 40' 30' 20' 

Right Ascension (J2000) Right Ascension (J2000) Right Ascension (J2000) 



3CJ.8'6 





7^145^^20' 10' 45'"0' 50' 40' 30' 20' gVg^'O' 30' 48"'0' 30' 47"'0' 30' 46"'0' 30' 1 o^24"'0' 50' 40' 30' 20' 



Right Ascension (J2000) Right Ascension (J2000) Right Ascension (J2000) 





h^gm^QS ^gniQS ^qS -^qS 

Right Ascension (J2000) 



66 02 

1 1^41 ""20' 41 ""0^ 40"'0^ 40' 20^ 

Right Ascension (J2000) 




11^55"'40' 30' 20' 10' SS'^O' 50 
Right Ascension (J2000) 



CLjT22^^332' 



12^27'"20' 10' 27"'0' 50' 40' 
Right Ascension (J2000) 





1 


40' 






' ^42' 














344' 


^^^^ 




1 






2 Q 48' 






-1 r5o' 




1 


- 



13"11'"20' 10' 11 ""0' 50' 
Right Ascension (J2000) 




13V7'"50' 40' 30' 20' 

Right Ascension (J2000) 
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Figure 1. Chandra images in the 0.7-7 keV energy range. The color bars reflect the 
number of counts detected by Chandra. SZ contour levels are (+2,-2,-4,-6,-8,...) times 
the rms noise in the short baseline data, after removal of radio sources; solid contours 
are for negative levels, and dashed contours are positive levels. The elliptical Gaussian 
approximation to the synthesised beam of the SZ observations are shown in the lower 
left corner. 
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3. Analysis of the SZA and Chandra data 



3.1. Models for the thermodynamic quantities 



We analyze the SZ and X-ray data using the Bulbul et al. (2010) model, which 



uses a consistent parameterization of the electron density temperature and pressure, 
related through the ideal gas law at all radii, i.e., Pe{r) = ne{r)kTe{r) for pressure p^^ 
electron density ng, and temperature Tg. All thermodynamic quantites depend on the 
gravitational potential, 

1 (l + r/r,)^-2-l 



(r) 



(3) 



_(/5-2) r/r,(l + r/r,)/5-2 
in which /3 describes the slope of the matter density at large radii and is a scale 



radius. The parameterization of the Bulbul et al. (2010) model does not allow the 



inner slope of the matter density to vary, which is fixed at r ^ as in the Navarro et al. 



(1997) model. The resolution of our SZ data can only effectively constrain the matter 



distribution on scales larger than the synthesized beam, which is of order 1 arcmin for 
these observations, and therefore we would not be able to place significant constraints 



on the inner slope. As explained in Bulbul et al. (2010), the potential is continuous at 
;5 = 2, the value of the Navarro et al. (1997) mass density model. The radial electron 



temperature profile is given by 

Te{r) = TQ(j){r)Tcooi{r), 



(4) 



where Tcooii"^) is the Vikhlinin et al. (2006) phenomenological core taper function. 



required to fit cool-core clusters, which is equal to one at large radii. The density 
is parameterized as 

ne(r) =neo0(r)V-^^(r). (5) 

in such a way that the pressure distribution is not altered by the presence of the cool 
core. At large radii, where the effect of the cool core vanishes, the thermodynamic 
quantities are related by a simple polytropic equation of state. The electron pressure 
profile is therefore parameterized as 

Pe{r) = Peo0(r)"+^ (6) 

and is independent of the presence of a cool core. The model therefore has five 
independent parameters for non-cool-core clusters: the scale radius r^; the index (5^ 
the polytropic index n, and the normalization constants for the three thermodynamic 
quantities which satisfy UeokTo = P^q. For cool-core clusters, the TcooI function 

+ {r/rcooiV 



adds three additional adjustable parameters. 



(7) 



To test for model-dependent biases, we also use the Arnaud et al. (2010) model 



to fit the SZ data. This model describes the cluster pressure profile using an analytic 
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function motivated by numerical simulations (Nagai et al. , 2007) and X-ray observations 
of the REXCESS sample, 

Pe,i 



Pe{r) 



{r/rpY [1 + {r/rpY 



\(b-c)/a ■ 



(8) 



The parameters p^^i and are left free in our fits to the SZ effect observations. The 
values (a, 6, c) are the power law indices that describe the (intermediate, outer, inner) 
slopes oipe{r). We use the "universal" values (a,6, c) = (1.05,5.49,0.31) obtained by 
Arnaud et al. ( |2010 ) from a fit to X-ray observations of the REXCESS sample. Note 



that Arnaud et al. (2010) find different best-fit values for cool-core clusters. We choose 



to use the parameters fit to the entire sample because our sample was not selected based 
on the presence of a cool core, and in fact contains a few non cool-core clusters, namely 
3C186, MS1137.5+6625 and CLJ1226.9+3332. 



3.2. Method of analysis 



As in previous work with the SZA (e.g., Mroczkowski et al. , 2009; Hasler et al. , 2011 ), we 



relate the point-source subtracted interferometric SZ visibilities to the unitless integrated 
Compton y by introducing Y{u^v)^ defined as 

V^{u,v) 



Y{u,v) = .... 

g{x) lo 

g{x) corrects for the frequency dependence of the SZ fiux, and /q 



(9) 



Here 

2{kBTcMBY / {hcY is the primary CMB intensity. The SZ models and compact radio 
sources are fit directly and simultaneously in Fourier space, where the statistical 
properties of the model fits are better understood and the noise is Gaussian. This 
is done simply by building up the sky brightness image, Fourier transforming it, and 
computing the likelihood of the model. 

The X-ray data consist of spectroscopic temperature measurements taken in cluster- 
centric annuli, and an X-ray image in units of surface brightness (counts s"-*^ cm~^ sr"-*^). 
The X-ray surface brightness Sx varies with the line of sight integral of the electron 
density and temperature distributions as 

1 



Sx 



niAee{Te,A)d£, 



(10) 



47r(l + z)^ 

where £ is the line of sight through the cluster, rie is the electron density, Tg is the electron 
temperature, A is the metallicity, and A^eiTe^A) is the X-ray cooling function (in units 
of counts cm^ s~^) as a function of electron temperature and metallicity. Each cluster 
was divided in a number of annuli according to the total number of photons detected, 
and for each annular region the temperature and abundance were free parameters. The 
surface brightness is only marginally sensitive to the choice of outer limit of integration 
in Equation 10: we find that the masses vary by less than 1% when the outer limit 



ranges between 2 and 5 Mpc. We therefore choose a limit of 2 Mpc, which corresponds 
to approximately the virial radius for clusters in this mass range. We use the |Mazzotta| 
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et al. (2004) definition of spectroscopic temperature in the comparison of model and 



observed temperatures in each annulus. 

We first estimate the pressure profile of the ICM by jointly fitting the SZ and X- 
ray data with the Bulbul et al. (2010) model. Both datasets are used simultaneously 
to constrain all three thermodynamic quantities, with the global shape parameters 
n and (and the cool-core parameters when applicable) linked among the profiles. 
Both datasets contribute to the determination of the shape of the pressure profile, 
with SZ observations contributing primarily at the largest radii where the sensitivity of 
Chandra to the diffuse cluster emission is limited. Instead of linking the normalization 
of the pressure profile (Peo) to the product of the normalizations of the density and 
temperature {n^o and Tq), we let the normalizations be free, and check a posteriori that 
Peo = Ueo X /cTq iu accordaucc with the ideal gas law. The normalization of the pressure 
is determined by the SZ data, and the normalizations of temperature and density by 
the X-ray data. 

This method results in the measurement of the shape of the pressure profile, 
Pe{r) / Peo, and two normalizations determined independently by each of the two datasets. 
The two normalizations are left free to vary because in principle systematic uncertainties 
in the two datasets could lead to different values, and we do not want to assume an a 
priori agreement between them. The fit uses a Markov chain Monte Carlo method 
(Bonamente et al., 2004), and computes the angular diameter distance assuming a 
Qa 0.73, Qm 0.27 and h 0.73 cosmology. 

To obtain a measurement of the integrated pressure that depends only on the SZ 



data, we also perform another fit in which we fix the shape parameters of the Bulbul 
et al. (2010) pressure profile to n = 3.5 and /3 = 2.0. These values correspond to 
the median of the values obtained from the joint fit. This pressure profile with fixed 



slope parameters is directly comparable to the universal pressure profile of [Arnaud et al. 



(2010), since both are determined by modelling of high-resolution X-ray data (from fits 
to the REFLEX sample for the [Arnaud et al.|[20T0| model, and from fits to the [Allen 



et al. 2008 observations for our model), and have just two free parameters (scale radius 



and normalization constant). In the following, we refer to this 2-parameter model as 
the Bulbul et al.[ ( [2010[ ) average pressure profile. 

Measurements of the ICM pressure using SZ and X-ray data are subject to different 
sources of systematic uncertainty that could affect the calculation of the l^ph parameter 
(Hasler et al., 2011). Systematic errors that integrate down with sample size include 
cluster asphericity, the effect of X-ray background, and the presence of kinetic SZ effect; 
these errors are included in the calculation of the ratio between the various measurements 
of ^sph (^500)5 and of the weighted averages and Xw values in Sections 4.1 and 4.2 
following the prescriptions of Hasler et al. (2011). 
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4. Integrated pressure measurements 



Joint SZ and X-ray fit using the Bulbul et al. (2010) model 



The integrated pressure, which we quantify in terms of the Compton y parameter, 
is expected to be a good proxy for total cluster mass. Since the SZA measures the 
integrated flux within Fourier modes on the sky, our SZ data relate most directly to the 
integrated Compton-y parameter l^yi. However, it is conventional in X-ray analyses to 
report spherically integrated quantities. We therefore quantify the integrated pressure 
using the spherically-integrated Compton y parameter l^ph out to r^oQ. The overdensity 
radius rsoo is given by 



rA = 



Mi 



tot\ 



(11) 



with A = 500, where Pc{z) is the critical density of the universe at the cluster redshift. 
The total cluster mass is calculated under the assumption of hydrostatic equilibrium; 
for the Bulbul et al.| (|2010) model, the total mass is given by 



MUr) 



1 



+ 



1/(1 - /5) - r/rs 



Tcoolir), 



(12) 



(/5-2) V/5-1 ■ {l + r/r^y-^ 
where the matter density normalization is given by pi = (/cTo(n+l)(/3— l))/(47rG/impr^); 
p is the mean molecular weight, and rUp is the proton mass. 



Using the method of analysis discussed in Section 3.2, we flrst compare l^ph 
normalized using n^o and Tq constrained by the X-ray data with l^ph normalized using 
Peo constrained by the SZ data. This comparison is summarized in Table [3| The 
normalizations are in good agreement: the weighted average of the ratio between the 
measurements using the SZ and X-ray normalization is 1.06 ± 0.04. This indicates 
that systematic uncertainties do not produce a large overall offset between the two 
observables. 

Below, we refer to l^ph as the measurement obtained from the joint flt using the X- 
ray normalization. We adopt this value since the joint proflle makes use of all information 
available from both the X-ray and SZ observations including the effect of the cool core, 
and since both normalizations are in agreement. 



4^2. SZ- only fit using the Bulbul et al (2010) average pressure profile 



We also flt only the SZA data to the Bulbul et al. (2010) average pressure proflle, which 



consists of the pressure proflle of Equation |6] with P^o and as free parameters and the 
two shape parameters flxed at n = 3.5 and /3 = 2.0. We use this model to compute Y 
as described above, which we refer to as l^ph,sz,Bio- The value of rsoo used in computing 
^ph,sz,Bio is determined from the joint flt. These results are shown in Table [4| and 
are plotted against the joint flt Fsph in Figure [2j We flnd that the weighted mean of 
the ratio between the measurements is given by l^sph,sz,Bio/^ph = 0.90 ± 0.05, where 
the uncertainty is the standard deviation of the weighted mean. A linear flt of the two 
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Table 3. Measurement of Integrated Fgph at rsoo from Joint X-ray and SZ Data Using 



the Polytropic Model 



Cluster 


^500 




-'^ spn 


(^500 ) 








SZ normalizat 


on 


X-ray normalization 


SZ-toX-ray ratio 




(") 


(10-11) 




(10-11) 




MACSJ0159.8-0849 


22i.i±l^:i 


8.30±g;^| 




9.67±i:it 


0.86±0:0f 


Abell 383 


268.5±i2;i 


4 Q9+0.77 




4.19±0:?§ 




MACS0329.7-0212 


138.4±if;7 


o no+0.53 




2 «l+0.59 




Abell 478 


7i4.3±it:i 


49.61±3;15 




60.62±|f4 


0-82±0:0f 


MACSJ0429.6-0253 


182.3±it:f 


9 7c:+0.49 
Z. ^O^Q 43 




3.30±0:ii 


0.83±0:it 


3C186 


72.i±|:| 


i.oi±g:ig 




0.86±oi| 




MACSJ0744.9+3927 


120.3±f:f 


5.04+0:66 




3.80±0:P 


i-33±g:lt 


MACSJ0947.2+7623 


196.2±ii:i 


1Q+0.73 




6m±\\l 


0.86±0:14 


Zwicky 3146 


265.7±i:^ 


12.17+1-22 




10.56±°:i 




MACSJ1115.8+0129 


200.0±^567 


7-74±0:i 






i-24±g:il 


MS1137.5+6625 


78.8±|-f 






0.73±0:ii 


1 4Q+0.26 


Abell 1413 


454.4±i0:3 


17 1 + 2.64 




23.12+2;32 




CLJ1226. 9+3332 


109.4±|:3 


o O1+0.34 
o.oizCq 34 




3.06±0:i| 




MACSJ1311. 0-0311 




2-36±0:ti 




2-3i±g:if 


1 00+0.26 


RXJ1347.5-1145 


218.0±i:6 


14.02±g:?i 




21.59±1J2 


n ^c^+0.04 


Abell 1835 


370.7±7.6 


31-41±Hi 




29.67±J-p 




MACSJ1423.8+2404 


i89.i±lt:t 


9 1 c+0.45 




2-52±g:iI 


0.86±0:?i 


MACSJ1427.3+4408 


i5o-5±4:i 


3.39±0:iS 






n 72+0.11 
' ^^0.11 


RXJ1504. 1-0248 


326.7±i2gi 


15.73±i:i 




18.03±i:5| 


n 07+0. 06 
^0.06 


MACSJ1532. 9+302 


i89.i±9:f 


c: n4_L0.65 




4-6i±g:iJ 




MACSJ1621. 6+3810 




9 c:o+0.29 




2-76±g:i 


0.93±0:12 


Abell 2204 


504.6±Jf:| 


44.97+2;99 




43.93±i:08 


i-02±g:gi 


MACSJ1720.3+3536 


170.5±|:f 


3.89+q"29 




3.93±0:t| 


0.98±0:io 


RXJ2129.6+0005 


297.6±ii:i 


10.78+i:g4 




10.48±i:i 


-1 04+0-10 
'■■"^^0.11 


Abell 2537 




7 27+0.81 

' -^'^0.77 




7-37±0:iI 


n qQ+0.09 







measurements to Sl y = x model results in a x 
we measure a scatter of 16%. 



2 

min 



35.3 for 25 degrees of freedom, and 



4^3. Comparison between the Bulbul et al. (2010) and Arnaud et al. (2010) pressure 
profiles applied to the SZ data 



The SZA data were also fit to the Arnaud et al. (2010) model using the same value of 



rsoo as above. The best-fit parameters are shown in Table [5| We compare the results 



from the Bulbul et al. (2010) average pressure model with the Arnaud et al. (2010) 



model in Figure [3j and find very good agreement: the weighted average of the ratio 



between the Bulbul et al. (2010) and the Arnaud et al. (2010) models is 1.05±0.06. 



A fit of the two measurements to a ^ = x model assuming the values are independent 
results in a Xw^^.G for 25 degrees of freedom, consistent with the presence of negligible 
scatter between the two measurements. The low value of Xw hkely due to correlated 
errors, since the two measurements make use of the same data. Figure [4] shows the 



average Arnaud et al. (2010) and Bulbul et al. (2010) pressure profiles for our sample. 



The two parameterizations result in fits that are consistent at all radii within r^oQ. The 
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Figure 2. Integrated pressure (l^sph,sz,Bio) from SZ data plotted against integrated 
pressure (Fgph) from the joint analysis, both measured out to the same value of r^oo- 
The dashed line is the curve y = x. 



consistency between the pressure profiles and the integrated Y{r ^qq) values measured 
from the two models indicate that the choice of parameterization for the gas pressure 
does not introduce a significant bias in the calculation of the integrated pressure within 

^500- 



5. Discussion 

The agreement we find between SZ and X-ray measurements of the l^sph(^5oo) parameter 
is consistent with a simple scenario in which the SZ decrement and the X-ray emission 
from massive relaxed clusters originate from the same highly-ionized thermal plasma, 
with only small contributions from other possible sources of emission. This result is in 
agreement with earlier ^ 30 GHz SZ studies using the Owens Valley Radio Observatory 
(OVRO) and the Berkeley Illinois Maryland Array (BIMA) millimeter arrays, in which 



the same value of the gas mass fraction was measured using SZ and X-ray data (LaRoque 



et al., 2006) Our results also support the findings by Melin et al. (2011) and Planck 



Collaboration et al. (2011b) of an overall agreement between the two measurements of 



the thermal pressure. 

We find scatter between the SZA and Chandra Fsph estimates at a level of 16%. 
A possible source of systematic error that could give rise to this scatter, and that is 
particularly relevant to our measurements out to rsoo, is elongation of the cluster along 
the line of sight. We use spherically symmetric models in the analysis; an intrinsically 
prolate cluster (elongated along the line of sight), when fit to a spherical model, will 
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Ysph,sz(W Arnaud et al. (2010) 



Figure 3. Integrated SZ pressure Fsph,sz calculated using the Bulbul et al. (2010) 



model {y axis) and the Arnaud et al. (2010) model {x axis), from a fit to the SZ 
data. The value of r^oo was determined by the joint modelling of the SZ and X-ray 
observations, and it is the same for both measurements. The dashed line is the curve 



y 



X. 



have its X-ray surface brightness — and therefore the corresponding l^ph parameter 



underestimated with respect to the corresponding SZ measurement (e.g., Cooray, 2000 



De Filippis et al. , 2005; Ameglio et al. , 2007). This is due to the quadratic dependence of 



the X-ray surface brightness profile on the density, as opposed to the linear dependence 
of the SZ effect. Our sample has just three clusters with a statistically significant 
deviation from the l^sph = ^ph,sz line, but in the direction of l^sph/^sph,sz > 1, and 
therefore consistent with oblateness (compression along the line of sight) rather than 
prolateness. The fact that the Allen et"aL] (2008) sample of relaxed clusters is X-ray 
selected may lead to including preferentially oblate clusters as their surface brightness 
will be boosted. An alternative interpretation for the presence of scatter between the 
SZA and Chandra estimates of Y is that some of these clusters are disturbed and have 



undergone a recent merger, as is almost certainly the case for RXJ1347. 5-1145 (Mason 



et al., 2010; Johnson et al., 2011). A merger would result in clumping of the gas, and 



therefore an overstimate of the gas mass and Y from X-ray measurements, as suggested 
by |Simionescu et al. (2011) to explain the observations of the Persues cluster. Clumping 
would not affect the SZ observations, because of the linear dependence of the signal on 
density. 

The fit of the SZ data to the universal pressure profile of Arnaud et al. (2010), and 
to the average pressure profile based on the Bulbul et al. (2010) parameterization of 
the pressure, are statistically acceptable for all clusters, with a similar foi* the two 
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r/r 



500 



Figure 4. Average pressure profiles from SZ fits to Bulbul et al. (2010) model (dark 



grey area, blue line), and to the Arnaud et al. (2010) model (light grey area, red line). 
The lines are the median of the 25 best-fit distributions, and the error bands are the 
68% confidence level. 



models. The agreement between Ygph at rsoo using the two models indicates that the 
integrated pressure is not highly sensitive to (reasonable) choices of parameterization. 

We have adopted throughout our analysis the value of rsoo determined from the 
joint SZ and X-ray observations. In the absence of X-ray information, one may instead 



adopt a fiducial value of the gas mass fraction fgas to determine rsoo (e.g., Joy et al. 



2001; Bonamente et al., 2008; Mroczkowski, 2011), or other means based on SZ-mass 



scaling relations. The additional assumptions required to estimate r^oo from SZ data 
only will likely contribute additional scatter to the l^sph — ^sph,sz relation, when r^oo used 
to measure l^ph,sz is estimated directly from the SZ data. 

6. Conclusions 



We have presented the joint analysis of SZA and Chandra observations of the [Allen 



et al. (2008) sample of massive and relaxed galaxy clusters. We have collected sensitive 
SZ data for all clusters at declination > —15° with no significant contamination from 
foreground or intrinsic radio sources, for a total of 25 clusters in the redshift range 
0.09 < z < 1.06. We also used the X-ray imaging and spectroscopic Chandra data 
that are available for all clusters, and performed a cluster-by-cluster comparison of 
the integrated pressure. The l^ph value estimated from the joint SZ and X-ray data, 
and from the SZ data alone, agree within a few percent at rsoo, indicating that the 
SZ and X-ray signal from massive relaxed clusters is consistent with a common thermal 
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Table 4. Best-fit parameters for the fit of the SZA data to the |Bulbul et al.| ( |2Q1Q] ) 

average pressure model, and integrated Y parameter out to rsoo- 



Cluster 


Peo 


Rs 


^sph,SZ,B10 


ratio l^sph,sz,Bio/^sph 




(10~^^ ergs cm~^) 


(arcsec) 


(10-11) 




MACSJ0159.8-0849 


27.50l^^gf 


40.99^^^5^^ 


9.07+i-^^ 


0.94+0.23 


Abell 383 


69.55 -|^g 22 


22.13+|:^J 


4.14+°-g 


0.99+0.29 


MACS0329.7-0212 


342.0011^-^° 


8.58i?-f2 


2.65+°-*? 


1.02+0.32 


Abell 478 


29.21+111 




70.981142-1? 


1.17+0.39 


MACSJ0429.6-0253 






2.86l°-53 


0.87+0.29 


3C186 


44.60+«|J 


13.28+^-^^ 


0.82l°;| 


0.95+0.32 


MACSJ0744.9+3927 


45.46lt?-3J 


2Q 91 + 13.04 


5.50 Q g4 


1.45+0.43 


MACSJ0947.2+7623 


31.44123-66 


34.31+^^4^ 


6.03+°-90 


1.01+0.31 


Zwicky 3146 


103.20+1^11 


25.36^1:^3 


11.89+i-« 


1.13+0.27 


MACSJ1115. 8+0129 




16.69+3-27 


7 99+° *'2 

58 


1.28+0.30 


MS1137.5+6625 


24.5619^57^ 


21.20+^-^9 


0-99io-l5 


1.36+0.37 


Abell 1413 


I2.53I1I5 


147.101^^-j^ 




1.95+0.81 


CLJ1226. 9+3332 


140.50+^^-?3 


129.00+^-J^ 


90+0.34 


1.07+0.29 


MACSJ1311. 0-0311 




54.26+^5-^2 


2 64+°-94 


1.14+0.47 


RXJ1347.5-1145 


296.7ll^|;58 


16 14+1-96 

±0.±^_2 20 


1 70+0.78 


0.63+0.13 


Abell 1835 


53.08l}|;« 


47.82liiof 


28.70+^-90 


0.97+0.22 


MACSJ1423. 8+2404 


■^6 CI +67.61 


24.80136-43 


4 44+3.12 

^•^^-1.42 


1.76+1.03 


MACSJ1427.3+4408 


90.97lg;« 


1 CI +8.13 
±0.0±_2 Qg 


2 26+°-^^ 

^•^"-0.44 


0.48+0.14 


RXJ1504. 1-0248 




zz. yo_g 25 


I2.9il?;g 


0.72+0.18 


MACSJ1532. 9+302 




Q7'+28.77 
^^•^'^-12.09 


c 79 + 1-83 
^- '^-1.10 


1.24+0.42 


MACSJ1621. 6+3810 




41 13+14.72 


9C'+0.55 


1.18+0.35 


Abell 2204 


OA -12+6.11 


89 47+16-45 

-11.79 


59.98+9-65 


1.37+0.33 


MACSJ1720.3+3536 


70 ^Q+50.98 
' ^•'-'^-29. 90 


18 54+5-49 

±0.0^_4 26 


op:+0.42 
<J-oC>_Q 29 


0.98+0.24 


RXJ2129.6+0005 




"^fi oc: + 28.27 
^"•^^-16.83 


1 2 ^7+3-37 
±Z.O^ _2 24 


1.21+0.38 


Abell 2537 




00.±0_-,^4 26 


8 74+1-99 
^•'^-1.43 


1.19+0.35 



origin. We therefore confirm the findings of Mehn et al. (2011 ) and Planck Collaboration 



et al. (2011b), and find no evidence for the presence of significant sources of systematic 



uncertainty in the measurements of the ICM pressure from SZ and X-ray observations 
of massive relaxed clusters. 



We also determined an average pressure profile based on the Bulbul et al. (2010) 
model, with shape parameters (n = 3.5 and /3 = 2.0) determined by a joint fit to Chandra 
X-ray data and our SZA observations of the Allen et al. (2008) sample of massive relaxed 
clusters. We have shown that measurements of the radial profile of the pressure out to 



rsoo, and of l^ph,sz at rsoo, agree between the Arnaud et al. (2010) and the Bulbul et al. 



(2010) average pressure profiles out to rsoo- Our conclusions indicate that both models 
are adequate for describing cluster radial pressure profiles and measuring the integrated 
thermal energy content in relaxed clusters. 
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Table 5. Best-fit parameters for the fit of the SZA data to the Arnaud et al. (2010) 
model, and integrated Y parameter out to rsoo- 
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